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We combine density functional theory, x-ray magnetic circular dichroism, multiplet calculations, 
and scanning tunneling spectroscopy to assess the magnetic properties of Fe atoms adsorbed on a 
thin layer of MgO(lOO) on Ag(lOO). Despite the strong axial field due to the O ligand, the weak 
cubic field induced by the four-fold coordination to Mg atoms entirely quenches the first order orbital 
moment. This is in marked contrast to Co, which has an out-of-plane orbital moment of L^ = ±3 
that is protected from mixing in a cubic ligand field. The spin-orbit interaction restores a large 
fraction of the Fe orbital moment leading a zero-field splitting of 14.0 ± 0.3 meV, the largest value 
reported for surface adsorbed Fe atoms. 


The magnetic properties of individual atoms adsorbed 
onto single crystal surfaces are remarkable. Their mag¬ 
netic anisotropy energies are several orders of magnitude 
higher than in bulk m- The main cause are the crystal 
field provided by the surface and the large orbital mag¬ 
netic moment preserved by the low coordinated adatoms. 
This moment interacts with the ligand field of the adsorp¬ 
tion site and is coupled to the spin moment by spin-orbit 
coupling (SOC). A system where the adatom exhibits the 
full gas-phase orbital moment is Go on a few monolay¬ 
ers (ML) of MgO(lOO) grown epitaxially on Ag(lOO) [5]. 
The binding of Co to the underlying O is strong while it 
is very weak to its next-nearest Mg neighbors. This cre¬ 
ates an almost perfect axial crystal field for the Co atom 
enabling the maximum possible magnetic anisotropy en¬ 
ergy for a 3c? atom. Moreover, the MgO layer decouples 
the atoms from the conduction electrons of the under¬ 
lying metal substrate reducing spin-flip scattering and 
increasing the magnetic relaxation times. 

To understand how individual Fe atoms react to this 
peculiar adsorption environment, we studied the mag¬ 
netic properties of Fe/MgO by combining density func¬ 
tional theory (DFT), x-ray magnetic circular dichro¬ 
ism (XMCD), multiplet calculations, and spin excitation 
spectroscopy (SES) with the scanning tunneling micro¬ 
scope (STM) [S]. We note that Fe ions in bulk MgO 
have long been considered as a model system for under¬ 
standing the interplay between crystal field and SOC ef¬ 
fects in 3d metals [13 E]. This interplay gives rise to a 
manyfold of low-lying excited states that strongly affects 
the magnitude of the orbital and spin moments and de¬ 
pends in a subtle way on local symmetry deviations from 
the ideal cubic environment. The discovery of high tun¬ 
nel magnetoresistance (TMR) [T2[ [13] and perpendicu¬ 
lar magnetic anisotropy in Fe/Mg(100) [TUlin] interfaces 


have strongly renewed the interest in that system due to 
technological applications in magnetic recording [MU]. 
Understanding the delicate balance among crystal field, 
SOC, and orbital hybridization between Fe and MgO is 
thus crucial, also to optimize the magnetic properties and 
downscaling of TMR devices. 

We find the lowest-lying zero-magnetic-field excitation 
at 14 meV. This energy is significantly larger than previ¬ 
ous reported values for Fe |2] , but also much smaller than 
the zero-field splitting of 58 meV reported for Co |5]- 
This difference between both elements is linked to the 
orbital symmetry of the magnetic atom. The low-lying 
orbital states of Fe = ±2 are mixed in the ligand 
field of the four next nearest-neighbor atoms (Mg). This 
fully quenches the orbital moment. In contrast, Co has 
Lz = ±3 and thereby is protected from such mixing. The 
spin-orbit interaction restores a significant fraction of the 
orbital moment for Fe, while it leaves the ground state 
orbital and spin moments largely unchanged for Co. The 
resulting second-order orbital moment observed for Fe is 
the prevalent situation when atomic spins are incorpo¬ 
rated in low-symmetry bonding geometries at surfaces, 
in bulk, and in molecules na^. 

Figure [^a) shows an STM image of two individual Fe 
atoms deposited at « 8 K on one ML MgO(lOO) grown on 
Ag(lOO) [3TH131I1S]' the Fe coverages of 0.01-0.03 ML 
(one ML is defined as one Fe atom per MgO(lOO) unit 
cell) used in the present study, we observe isolated Fe 
atoms rather than clusters. We find only one Fe species 
with an apparent height of 180 pm. The lowest energy 
adsorption site found in DFT is the one on-top of oxy¬ 
gen 1^1211 US]- Figure [^b) reveals that the O beneath 
the Fe is displaced by 40 pm upwards from the MgO 
plane [21]. The ground state electron density shows a 
significant charge transfer between the Fe and the near- 
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Figure 1. (color online), (a) STM image of two Fe atoms 
on a ML MgO(lOO) grown on Ag(lOO) (4 nm x 4 nm, tunnel 
current 7t = 5 pA, tunnel voltage Vt — 100 mV), (b) Side 
view of DFT calculated binding geometry and charge density 
(color scale le/(au)®, Fe green, O red, Mg blue). Middle 
sketch shows top view ball model of the binding geometry, (c) 
Oblique view of DFT-calculated valence electron spin density 
contours (positive spin polarization - red, negative - blue). 


More insight into the magnetic levels and the evolu¬ 
tion of their energies is gained from multiplet theory |43) . 
The calculated spectra shown in Fig. [^are in very good 
agreement with experiment for both incident beam direc¬ 
tions. In addition, the experimental magnetization curve 
is perfectly reproduced by the line showing the out-of- 
plane projected field-dependent total magnetic moment 
{2Sz{B)) + {Lz{B)) derived from multiplet calculations. 
In these calculations we included charge transfer to the O 
ligand, leading to configuration mixing, the axial ligand 
field due to the nearest-neighbor O atom (Ds), the cu¬ 
bic one due to the four next-nearest neighbor Mg atoms 
(Dq), SOC (C), and the external magnetic field (B). Best 
agreement is obtained with a 90 % d 6 -b 10 % d7 ( con¬ 
figuration of the Fe atom, where I refers to a ligand hole 
in the neighboring O atom. 

The configuration mixing results in a 10-fold degener¬ 
ate ground state. Ds moves this state down in energy, 
while all other states move up. The resulting ground 
state ((Tz) = ± 2 ) 0 ((S'z) = ±1.96, ±0.98, and 0) is 
shown on the left hand side of Fig. that illustrates 
its evolution under the action of Dq, and B. In 
marked contrast to Co, Dq strongly perturbs the low¬ 


est O atom with an overall positive charge of ±0.44 e on 
the Fe [2^. Note that Co on the same binding site is 
nearly charge neutral [ 8 ]. 

With the generalized gradient approximation (GCA) 
and on-site Coulomb interactions {U = 3.2 eV) [37] for 
the Fe d-states, we find a total spin moment of 3.6 fiB on 
the Fe atom. The majority spin density, red in Fig. [^c), 
is mostly axially symmetric and exhibits an induced po¬ 
larization of the O atom slightly increasing the overall 
spin moment to 3.7 fj,B- The four-fold symmetry of the 
binding site shows up strongly in the minority spin den¬ 
sity (blue) in contrast to Co, which exhibits nearly per¬ 
fect axial spin density | 8 ]. The calculated spin of the 
Fe atom is to very good approximation independent of 
the MgO thickness. This facilitates the interpretation 
of synchrotron measurements on samples having several 
coexisting MgO thicknesses |25j . 

X-ray absorption spectra (XAS) and the resulting 
XMCD signal [55j are shown in Figs, [^a) and (b) (see 
Ref. for experimental details). The Fe ±3 and ±2 
edges exhibit sharp multiplet structure characteristic of 
an ensemble of single adatoms on identical adsorption 
sites [39]. The XMCD signal is larger at normal than 
at grazing incidence, revealing out-of-plane easy magne¬ 
tization axis. The sum rules |3DM3] yield a large orbital 
moment of (Lz) = 1-74 ± 0.11/ie and an effective spin 
moment of {2Sz) + {7Tz) = 2.46 ± 0.11 (assuming 
3.9 d holes, as found in the multiplet calculations below; 
Tz is the out-of-plane projection of the atomic magnetic 
dipole moment). The measured out-of-plane magnetiza¬ 
tion curve is shown in Fig. [^c). 
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Figure 2. (color online), (a) Measured and simulated XAS 
over the Fe L 3 and L 2 edges for 0.03 ML Fe on an MgO film 
on Ag(lOO) with an average thickness of 3 ML. Red curves 
are for normal (0 = 0 °) and blue for grazing (0 = 60°) inci¬ 
dence of x-ray beam and B held (T = 2.5 K, B = 6.8 T, total 
electron yield mode), (b) XMCD spectra for both geometries, 
(c) Out-of-plane magnetization curve measured by hrst sat¬ 
urating the sample at 6.8 T (red) and — 6.8 T (green) and 
then moving to the respective held value (T = 2.5 K). The 
solid line represents {2Sz{B)} + {Lz{B)) determined by the 
multiplet ht with a saturation moment of 5.2 fiB- (d) Sketch 
of the measurement geometry. 
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Figure 3. (color online). Energy level diagram resulting from 
the multiplet simulation of the XAS and XMCD spectra (for 
full diagram see [5^). The cubic crystal field quenches the or¬ 
bital moment and creates two spin quintuplets with Bi (blue) 
and B 2 (red) symmetry. SOC creates an essentially 2-fold de¬ 
generate ground state with largely restored orbital moment 
of (Lz) = ±1.25. The Sz and Lz values in the labels are the 
respective expectation values, the () signs have been omitted 
for brevity. The Zeeman energy splits these states into the 
five states |0) - |4). The consecutive states deriving from the 
B 2 manifold are labelled |5) - |9). 


(a) (b) 




Figure 4. (color online), (a) STM differential conductance 
(dl/dV) spectrum on Fe/1 ML Mg0(100)/Ag(100), dl/dV 
on bare MgO is shown for comparison (T = 0.6 K, B = 
0 T, K modulation with rms amplitude V)nod = 150 /rV and 
frequency / = 806 Hz, set point before opening the feedback 
loop It = 1 nA, Lt = 30 mV), (b) Positive conductance step 
at out-of-plane fields of 0 T (green) and 4 T (blue). Inset: 
Field splitting of conductance step energies, (c) SES feature 
corresponding to the superposition of V35, V25, V47, and V40, 
measured with a spin-polarized tip (T = 1.2 K, B = 2 T, 
Hmod = 1.5 mV, It = 1 nA, V) = 100). Inset: spin-polarized 
spectrum in the same energy window as (a), (d) Sketch of the 
magnetic states and the allowed SES excitations. 


est multiplet and creates two spin quintuplets with fully 
quenched orbital moments, ((Lz) = 0 ) 0 ((*S'z) = ±1.96, 
±0.98, and 0). The one with Bi symmetry is the ground 
state and drawn in blue. SOC splits it into three energy 
levels, where the lowest restores more than half of the 
free-atom orbital moment by coupling the two lowest or¬ 
bital levels in a second-order perturbation m- At zero 
magnetic field, the new ground states {Lz) = ±1.25 and 
(Sz) = ±1.96 are to a very good approximation two-fold 
degenerate (ignoring energy differences of a few neV). 
The magnetic field lifts the remaining degeneracy of the 
five states labelled |0) - |4). The excited spin quintu¬ 
plet (red) has B 2 symmetry, lies « 100 meV higher in 
energy, and has its orbital and spin magnetic moments 
anti-aligned. This results in smaller total magnetic mo¬ 
ments and hence a smaller Zeeman splitting. 

Our STM-SES measurements on individual Fe atoms 
determine the energy splitting of the lowest lying mag¬ 
netic states with high precision, and thereby complement 
XMCD. Figure |^a) displays clear conductance steps 
with 15 % amplitude located at ±14.0 ± 0.3 mV (the 
error bar refers to variations between atoms at differ¬ 
ent locations of the MgO film). The magnetic nature of 


the underlying inelastic excitations is demonstrated by 
the splitting of the excitation energy in an out-of-plane 
magnetic field shown in Fig. Bb). For in-plane fields 
the splitting is absent indicative of an out-of-plane 
easy axis. The observed zero-field splitting of 14 meV is 
more than twice the largest values seen for individual Fe 
atoms adsorbed on other surfaces [5] and our magnetic 
anisotropy energy [33] approaches that reported for Fe 
atoms in linear molecules [3H] . 

Connection between the transitions excited in SES and 
the states derived from the multiplet calculations is es¬ 
tablished by the level diagram shown in Fig. B^)- The 
horizontal position of the states |0) - |9) indicates their 
out-of-plane projected magnetic moments. The blue ar¬ 
rows are the first two excitations possible for tunnel elec¬ 
trons, |0) —)■ |2) and |1) —>■ |3), with the corresponding 
tunnel voltages labeled V 02 and V 13 . These are the exci¬ 
tations giving rise to the conductance steps in Figs.Qa) 
and (b). At zero field, V 02 = V 13 = 14.0 mV. This is in 
excellent agreement with the level separation of 13 meV 
found in the multiplet calculations. 

As expected from the Zeeman splitting in Fig. V 02 
and V 13 shift up and down symmetrically in an external 
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out-of-plane field, see Fig. [^b). As seen from the inset, 
the splitting is linear in field and amounts to 1.9±0.3 mV 
at 6 T. The amplitude of V 13 is significantly smaller than 
that of V 02 because the corresponding transition starts 
from the excited state |1). That V 13 is visible at all im¬ 
plies that the tunnel current drives the occupation of the 
| 1 ) state, and that its lifetime is longer than the mean 
time between tunneling electrons of the order of 1 ns [46] . 
From the slopes in the inset of Fig. [^b) we derive an ef¬ 
fective 5 *-value of 2.57 ± 0.06 in the z-direction. This is 
significantly above the free-electron value of = 2 and 
thus confirms the presence of a large unquenched orbital 
magnetic moment for Fe on MgO [13148] . From XMCD 
and the multiplet calculations we derive g* as the dif¬ 
ference of the out-of-plane projected magnetic moments 
between states |0) and | 2 ), g* = A{Lz) + 2A{Sz) = 2.46, 
again in very good agreement with SES. 

The magnitude of the zero-field splitting and of the 
orbital moment strongly depend on the cubic ligand field. 
However, the simulated spectra are not very sensitive to 
small changes of Dq. This illustrates the strength in 
combining SES and XMCD. The first measures the zero- 
field splitting directly and very precisely, while the second 
identifies the ground and excited states with their spin 
and orbital magnetic moments, as well as their angular 
anisotropies. 

The strong sensitivity to the cubic ligand field shows 
that Fe exhibits a fundamentally different magnetic be¬ 
havior on the MgO surface compared to Co on the same 
substrate | 8 ]. The key difference is that the cubic field is 
able to quench the Fe orbital moment because the four¬ 
fold symmetry can mix the = ±2 components, result¬ 
ing in (Lz) = 0 after the ligand field contribution. The 
Lz = ±3 components of Co are protected from such mix¬ 
ing, leading to a linear response to SOC | 8 ] . Fe represents 
the more common case [49], where SOC restores some of 
the otherwise quenched orbital moment, due to mixing 
with the next-higher orbital state. 

Spin-polarized (SP) STM tips reveal an additional con¬ 
ductance step at ^ 105 mV, see Fig. Qc). This step is 
absent for non-polarized tips. Therefore, we assign it 
to an electron driven occupation change of the magnetic 
states at the respective threshold energy (spin pumping). 
Evidently, this part of the conductance change in spin 
polarized spectra depends on the set-point current, since 
the tunneling electrons must arrive frequently enough to 
probe the excited states before they decay [46]. At the 
V 02 and V 13 step energy, inelastic spin excitations and 
spin pumping are superimposed creating the jagged edges 
shown in the inset of Fig. Qc). 

In order to identify which states contribute to the 
105 mV step, we note that the transition rate from the 
initial to a final state follows Fermi’s golden rule and 
is proportional to the interaction potential of the atom’s 
spin with the spin of the tunneling electron niisiisniisii- 
This model implicitly requires that the orbital moments 


remain constant, while the tunneling electrons can ex¬ 
change spin angular momentum with the atom according 
to the usual selection rule ASz = ±1, 0. Note that the 
selection rule of the orbital moments is fulfilled if part of 
the initial orbital wave function overlaps with the final 
orbital wave function. This is the case for all lowest 10 
states of Fe, as they all have components of Lz = ±2, de¬ 
spite the changing expectation value (Lz). Using these 
selection rules, we can assign the ^ 105 mV step to the 
four transitions V 35 , U 25 , V 46 , and V 47 that are very close 
in energy [52] . This again implies lifetimes r > 1 ns of 
one or more of the initial states |2), |3), and |4). The mul¬ 
tiplet calculations estimate V 35 ~ 134 meV and predict 
a transition probability one order of magnitude smaller 
than for the V 02 transition [^. The small intensity ex¬ 
plains the absence of this step in non-polarized dl/dV 
spectra. The energy is in reasonable agreement with the 
one directly probed by STM. 

Transitions between states belonging to different spin 
multiplets have been observed in spin chains {e.g. sin¬ 
glet to doublet and triplet states) [33 [i3| and molecular 
magnets m- Unique to Fe on MgO, the excitation from 
lower to upper spin multiplet involves a transition from 
aligned (lower multiplet) to anti-aligned (upper multi¬ 
plet) orbital and spin moments. 

The present study reveals how the perpendicular mag¬ 
netic anisotropy and enhanced orbital moment of Fe 
atoms on MgO arise from the interplay of SOC and the 
low-symmetry ligand field at O sites. Despite the simi¬ 
lar chemical environment, we find important differences 
between Fe and Co adsorbed onto MgO thin films. For 
both elements, the axial ligand field of the oxygen fully 
preserves the gas phase orbital moments of L = 2 , respec¬ 
tively, L = 3. However, the cubic ligand field created 
by the much weaker interaction with the next-nearest 
neighbors (Mg) acts very differently on both elements. 
For Fe it quenches the orbital moment due to mixing 
of the Lz = ±2 components of the orbital wave func¬ 
tion, while the Lz = ±3 components of Co are protected 
from that mixing by symmetry. In the case of Fe, SOC 
restores a very large fraction of the orbital moment of 
(Lz) = 1.25/tb, leading to a high zero-field splitting 
of 14 meV unprecedented for Fe adatoms. The mag¬ 
netic levels derived from XMCD and multiplet calcula¬ 
tions agree extremely well with the excitation energies, 
their field splitting, and excitation amplitudes observed 
in STM-SES. Remarkably, spin-polarized STM can also 
probe excitations to higher spin multiplets. 
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